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We study the dynamics of SU(2) gauge theory with Nf = 6 Dirac fermions by means of lattice 
simulation to investigate if they are appropriate to realization of electroweak symmetry break- 
ing. The discrete analogue of beta function for the running coupling constant defined under 
the Schrodinger functional boundary condition are computed on the lattices up to linear size of 
L/a — 24 and preclude the existence of infrared fixed point below g^ ^ 7.6. Gluonic observables 
such as heavy quark potential, string tension, Polyakov loop suggest that the target system is in 
the confining phase even in the massless quark limit. 
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After the seminal work [|l]], great attention has been paid to the possibility to study nearly 
conformal dynamics of gauge theory by means of lattice simulation, which are expected to trigger 
electroweak symmetry breaking. So far, many works have been involved in the calculations on 
SU(3)c gauge theory with Np Dirac fermions in the fundamental representation. Recently, it was 
pointed out that three-color QCD with ten-flavors may be conformal in the infrared limit with 
large mass anomalous dimension [|2[ ^. However, large volume simulations indicated that chiral 
symmetry breaking occurs in the twelve-flavor system The large scale simulation via Monte- 
Carlo renormalization group [Q] will be effective to extract the genuine dynamics of the system. 

Here we focus on a series of SU(2)c gauge theories with Np Dirac fermions in the fundamental 
representation, and try to find out a gauge system with the critical number of flavors A^"', at which 
the chiral condensate "i/7i/A gets large anomalous dimension. SU(2) has a symplectic form so that its 
fundamental representation is pseudo-real, and that the chiral symmetry of A^f -flavor system is en- 
hanced to S\]{2Nf) D SU(A^f )l X S\]{Nf)r x U(1)v- If chiral symmetry is spontaneously broken, 
the plausible unbroken symmetry subgroup is SP(2A^/r). But, this enhanced unbroken symmetry 
may contain the electroweak symmetry Gw = SU(2)lxU(1)y depending on the representation 
of new "quarks" under Gw- It is thus inevitable to examine the vacuum alignment issue [^, i.e., 
whether Gw is broken or not. The transition between confinement and deconfinement is argued 
to be second-order in pure Yang-Mills theory with SU(2)c while it is first-order for SU(3)c- 
Therefore, SU(2)c chiral dynamics can differ even qualitatively from those of SU(3)c, in particular 
at A^^". 

Actually, two-color QCD has also been studied thus far. Iwasaki et al. showed that A^^ = 
3 system is conformal in the infrared (IR) limit ^ through study of phase structure of Wilson 
fermions, while the perturbatively calculated j3 function suggests that 6 < A^^"^ < 8. Afterwards, 
running gauge coupling constant has been calculated nonperturbatively for two-color QCD with 
six-flavors [^-11], and eight-flavors [p^, implying that those systems are conformal in the IR 

limit ig In], |1|]. 

We note that SU(2)l gauge theory with three generations of quarks and leptons is exactly the 
system of our concern here, the two-color QCD with six massless Dirac fermions. Our question 
is if the quantum-mechanical dynamics of the presumed copy of the existing gauge symmetry and 
fermionic matters but with large A 'qcd" could play the role of triggering spontaneous breakdown 
of the electroweak symmetry, SU(2)l x U(1)y- 

The purpose of this article is to report our first result for the dynamical features of two- 
color QCD with six-flavors according to the simulation in the framework of lattice gauge the- 
ory. Throughout this work, the standard Wilson plaquette gauge action with unimproved Wilson 
fermions is used for simulation. 

We measure gsF{L/a, g^) for sets of {L/a, = ^) defined under the Schrodinger functional 



boundary condition with the twist angle for the quark fields set to [ 1 3 ] . Compared to the preceding 
work [^, computation on larger lattices, L/a = 6, 8, 12, 16, 18 and 24, is done with fine tuning 
of the critical value of hopping parameter for every pair of {L/a, g^). Data are fit to the presumed 
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functional form 
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where the coefficient Pija.i is obtained by the one-loop calculation. Using the fit result gp)' 
the discrete beta function [fi^ 
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for a fixed step scaling factor s can be calculated. 

Figure |l| shows our result for its contin- 
uum limit B{u, s = 2) taken in the same man- 
ner as in Ref. [[l5|]. The result precludes the 
existence of infrared fixed point below 7.6. 
We are also trying to calculate the anomalous 
dimension of the chiral condensate YY- We 
find that the result is roughly consistent with 
the perturbative prediction and the systematic 
uncertainties are under examination. 

The above result suggests that A^f = 6 is 
not far from A^"' and motivates us to investi- 
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Figure 1: Discrete beta function for the step scaling 
s — 2. The horizontal axis is the inverse of squared 
renormalized coupling constant. Red plots were ob- 
tained from all data, while blue plots were obtained 
from data without L/a = 24. 



gate the theory further from a different point 
of view. We thus study the dependence of 
mesonic spectrum on quark masses. First, in 
order to see the phase structure as a statistical 
system and fix the simulation parameters, we performed a scan on the (j8, K")-plane with relatively 
small lattices. Figure ^ shows the plaquette (W) as a function of I/k" on 8^ x 24 and 8^ x 32. For 
j8 < 1.7, there is a region of K in which (W) changes rapidly. 

Thus, we choose j8 = 2.0 to simulate larger lattices for our target system. 

Figure ^ shows the masses of the lightest pseudoscalar and vector mesons, where the depen- 
dence on K is translated to that on ampcAC via Figure ^ and the following properties are observed: 

(1) The meson masses are bounded from below at small ampcAC- The lower bound anf^^{L/a) 
depends on L/a. 

(2) Even for small ampcAC, the ratio of the mass of the vector meson to that of the pseudoscalar 
meson is not far from 1 . Figure || shows that the scalar meson is paired with the axial-vector 
meson ' , and that their masses are well above those of the lightest pseudoscalar and vector 
mesons before saturation occurs. 



' The fit is performed for the two-point correlation functions in the scalar and axial- vector channels before substantial 
fluctuation sets in at large t/a. 
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Figure 2: Plaquette, lyV), versus 1/k" on the lattices 
of spatial linear size L/a — % with )3 = 1.5, 1.6, 1.7 
and 1.9. 
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Figure 3: Lightest pseudoscalar (o) and vector (•) 
meson masses versus PCAC mass mpcAC in lattice 
unit at /3 = 2.0 on four different volumes. 
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Figure 4: PCAC mass mpcAC versus 1/k at j3 = 
2.0. The inset zooms in the vanishing mpcAC region. 
Finite size effect is negligibly small. 
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Figure 5: Masses of lightest pseudoscalar (P), vec- 
tor (V), scalar (S) and axial-vector (A) mesons for 
L/a = 24 and L/a = 32. 



For L/a = 8, we also checked that the property (|T|) persists for Iwasaki gauge action and/or clover 
fermions, and that am^'(L/a = 8) depends modestly on the types of action and j3. This saturation 



phenomenon was also observed in SU(3)c gauge theory with two sextet quarks [|16|]. When we 
look at Figure ^ from large ampcAC, the meson masses branch at larger ampcAC with am^^{L/a) 
for smaller L/a, which strongly indicates that this saturation originates from finite size effect. 

We comment on the dependence of anf^^{L/a) on the linear size L/a of the system. The spatial 
correlation length cannot become larger than the system size L/a and saturated at L/a. Thus, 
if the Compton wavelength 27r/(am,^) coincides with ^^f,, ai7f^^{L/a) then decreases with L/a, 
accounting for the dependence of anf^^{L/a) on L/a in Figure However, this is actually not 
the case in SU(2)c gauge theory with two adjoint Dirac fermions, where amp\L/a) increases with 
L/a [p^. We recall that the finite size correction consists of two terms [|I9|]; the term induced 
by PP scattering which increases the pseudoscalar meson mass amp, and the one induced by the 
propagation of 0++ which decreases amp. Since the two adjoint fermion system contains a glueball- 



^The finite size effect observed here differs from the power correction found in Ref. [p^. 
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rich 0++ lighter than P [20], the latter contribution is possibly larger than the former. The decrease 
of anif^{L/a) with L/a in our target system suggests that the PP scattering contribution is the 
dominant source of the finite size effect. The knowledge on the mass of the lightest 0++ will surely 
give us more comprehensive understanding on this issue. 

Our interest is if the property reflects the dynamics of two-color QCD with six flavors of 
almost massless quarks. The high degeneracy in masses between vector and pseudoscalar mesons 
was also observed in the SU(2)c gauge theory with two adjoint fermions [|l8]|, which is considered 
to be conformal in the deep IR. Meanwhile, such a degeneracy reminds us the spectrum of bound 
states of massive quarks, lUq » Aqcd- Another possibility is that it may occur as a consequence of 
finite size effect; even if the theory is confining, the system size is too small for the confining force 
to act between quark and anti-quarks so that they are bounded solely by Z-copies of Coulombic 
forces. The knowledge on the heavy quark potential will help to answer these questions. 
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Figure 6: Heavy quark potential obtained from the lattice with 24^ x 48, j3 = 2.0, 0.1490 < if < 0.1498 (left 
panel), and 0.1498 < K" < 0.1502 (right panel). Note that the scale of the ordinate of the right differs from 
that in the left. 
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The heavy quark potential is ex- 
tracted from the measurement of various 
sizes of Wilson loops. To reduce short dis- 
tance fluctuation, four-level Wilson flow 
smearing |21] with the weight e = 0.01 in 
the exponent was performed for link vari- 
ables in the spatial directions. Figure ^ 
plots the heavy quark potentials on the lat- 
tices with L/a = 24 and various k. See 
Figure ^ for the correspondence of 1 / K" 
and ampcAC- We can see that the heavy 
quark potential contains the component of 
linear term, i.e., confinement. 

Figure |7| shows the squared root a^/a 
of the string tension. The comparison of a^/a for L/a = 24 with that for L/c? = 32 indicates that 
the finite size effect is not substantial for string tension. Moreover, there is no such a supporting 
evidence that it approaches to zero in the limit ampcAC — ^ 0, suggesting that the theory is confining 
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Figure 7: Masses of lightest pseudoscalar (P) and 
squared root of string tension a^/a. 
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in the chiral limit. However, am^ relative to a^/o in Figure ^ cautions that if the system exhibits 
spontaneous chiral symmetry breakdown, the simulation on the lattice, say with size L/a = 48, will 
be necessary to capture genuine chiral dynamics. 
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Figure 8: Distribution of Polyakov loops in the Figure 9: Distribution of Polyakov loops in the di- 
z-direction {Nz = 32, periodic) on the lattice with rection (Nt = 8 and anti-periodic) on the lattice with 
32^ x64, j3 =2.0, fc^ 0.14965. 18^x48 x 8,j3 =2.0, fc^ 0.1492. 



Lastly, we study the distribution of Polyakov loop P in the thermal direction. P takes its value 
in a real number in the SU(2) gauge theory. Its vacuum expectation value (P) is not a good mea- 
sure of confinement/deconfinement transition in presence of dynamical quarks in the fundamental 
representation, which explicitly breaks the center Z2. In spite of this fact, its distribution will help 
to capture the essence of the dynamics; the system is deconfining if the distribution of P does not 
cover and is completely asymmetric with respect to the origin, while it is confining if it peaks at 
the origin and distributes almost symmetrically. 

We can see that the periodic boundary condition is equivalent to the anti-periodic boundary 
condition in SU(2)c gauge theory if all matter fields belong to the representations of Z2-odd con- 
jugacy classes. (Only the sign of a Z2-odd observable depends on the condition.) Thus, we regard 
one of the spatial directions with 32 sites on lattices of 32^ x 64 as a thermal direction and measure 
P along this direction. Figure |8| shows the distribution of P for K" = 0. 14965. The periodic/anti- 
periodic boundary condition will reduce the dynamical modes so that the finite size effect will act 
to order the system ((P) 7^ 0, i.e., tendency of deconfinement). Even though the aspect ratio is 
32/32, the distribution of P in Figure || exhibits disordering, suggesting confinement. In contrast, 
the distribution of P for A^r = 8 in Figure ^ does not cover 0, suggesting deconfinement. The 
question whether this transition is really thermal type needs further study. 

To summarize, we show some evidence supporting confinement in the two-color QCD with 
six-flavors of quarks. As can be seen in the inset of Figure ^, the mass splitting between pseu- 
doscalar meson and vector meson increases gradually in the small quark mass region until satura- 
tion is encountered. The study with larger size of lattices is necessary to approach smaller nip / ^fc 
and to get a definite conclusion. Now, it is interesting to start the calculation of the other physical 
quantities, such as the mass of the lightest particle in the 0++ channel relative to the decay constant. 
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